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Abstract. Consider the semiclassical limit, as the Planck constant h — > 0, of 
bound states of a quantum particle in a one-dimensional potential well. We 
justify the semiclassical asymptotics of cigcnfunctions and recover the Bohr- 
Sommerfeld quantization condition. 



1. Introduction 

1.1. Wc study the limit as h — » of eigenfunctions ip(x) = ip(x;\,h) of the 
Schrodinger equation 



-h 2 ip"(x)+v(x)ip(x) = \ip(x), v(x) = v(x), ipeL 2 (R), (1.1) 

for A close to some non-critical energy Ao (that is v'{x) ^ for x such that v(x) — 
Ao). We assume that the equation v(x) = A has exactly two solutions (the turning 
points) x± = x±(X) and that v[x) < A for x G (x-,x+). Thus, (x-,x+) is a 
potential well and the energy A is separated from its bottom. We suppose that 
eigenfunctions ip(x) are real and normalized, that is 

ip 2 (x)dx = 1. 

It is a common wisdom that the limit of ip( x ) = VK 3 -! \ ^) as ft — >■ is described 
by the Green-Liouville approximation away from the turning points x±. In neigh- 
borhoods of the turning points the asymptotics of ij)(x) is more complicated and 
is given in terms of an Airy function. Surprisingly, we have not found a precise 
formulation and a proof of this result in the literature. Our goal is to fill in this 
gap. We follow here the scheme suggested by R. E. Langer and thoroughly exposed 
by F. W. Olver in his book [7J. 

The detailed asymptotics of ^(x) described in Theorems 12.51 and 14.41 allows one 
to recover the classical Bohr-Sommerfeld quantization condition on A (see Theo- 
rem [47TJ) . Actually, we prove somewhat more establishing a one-to-one correspon- 
dence between eigenvalues of the Schrodinger operator H% = —h 2 d 2 /dx 2 + v(x) 
from a neighborhood of a non-critical energy and points (n + 1/2)H where n is an 
integer. This implies the semiclassical Weyl formula for the distribution of eigen- 
values of the operator Hn as H — ¥ with a strong estimate of the remainder. It 
turns out (see Corollary I4.2[) that this remainder never exceeds 1. We also obtain 
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in Theorem 15.21 the quantization condition for discontinuous functions v(x). This 
formula generalizes that of Bohr and Sommerfeld and is probably new. 

We note that the Bohr-Sommerfeld quantization condition is also well known in 
much more difficult multidimensional problems. In this context we mention book 
P] by V. M. Babich and V. S. Buldyrev (where the ray approximation is used), book 
[3] by M. V. Fedoryuk and V. P. Maslov (where the Maslov canonical operator is 
used) as well as papers [5] by B. Helffer et D. Robert and [4] by B. Helffer, A. 
Martinez and D. Robert (where the methods of microlocal analysis are used). 

However, in the one-dimensional problem it is more natural to rely on methods of 
ordinary differential equations. Such an approach was developed by M. V. Fedoryuk 
(see his book [2]) for analytic potentials. In this case one can avoid a study of 
turning points so that the Airy function does not appear. 

1.2. The asymptotics of eigenfunctions yields (see Proposition l4.5|) asymptotics 
of observables 

f oo 

w(x)ijj 2 (x;X,h)dx (1.2) 



for sufficiently arbitrary functions w(x). For example, we can take for w(x) char- 
acteristic functions of Borel subsets of R or choose w(x) — v(x). This gives the 
asymptotics of the kinetic energy 

/>oo 

K{\h):=h 2 tl>'(x;\,h) 2 dx = K c i(\)+0(h 1/3 ) (1.3) 



as h — > uniformly for A from a neighborhood of the point Ao. The leading term 
K c i(X) (the index "d" stands of course for the corresponding classical object) is 
given by the expression 

K cl {\)= {X-v(x))^ 2 dx( (A - v(x))~ x / 2 dx) . (1.4) 

Ji_(A) X Jx-{\) J 

Note that the integrals here are taken over the classically allowed region and that 
K c i(X) coincides (see subsection 4.3) with the averaged value of the kinetic energy 
of a particle of energy A in classical mechanics. 

We emphasize that our derivation of the quantization condition and of asymp- 
totic formulas for observables (|1.2j) requires Airy functions although they do not 
enter into the final answer. However, we do not know how to avoid Airy functions 
without additional assumptions on v(x). 

2. Semiclassical solutions of the Schrodinger equation 

2.1. It is convenient to rewrite equation (|1.1|) as 

- <(x) + tT 2 q{x)u h {x) = 0, (2.1) 

where 

q(x) — q(x; A) = v(x) — A. 

We need some regularity of the function v(x) and a weak condition on its behavior 
at infinity. 

Assumption 2.1. The function v 6 C 2 (R) and, for some po > 1, the function 

(\q(x)\- S q'(x) 2 + q(x)- 2 \q"(x)\) f \q(y)\^ 2 dy *' 

Jo 
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is bounded for sufficiently large \x\. 

The last condition is satisfied in all reasonable cases. For example, if v(x) — > 
vq > A, it is sufficient to require that 

v'(x) 2 + \v"(x)\=O{\x\- p0 ) : p >l, |a;|-). oo. 

It is also satisfied if v(x) behaves at infinity as |x| Q or e a ^ where a > 0; in these 
cases po = 2. 

We consider the case of one potential well. To be more precise, we make the 
following 

Assumption 2.2. The equation v(x) = A has two solutions x + = x+(A) and 
x_ = X-(X). We suppose that X- < x+, v{x) < A for x € (x_,x+), v(x) > A for 
x $ [x_, x+] and 

liminf v{x) > A. 

|a;| — >oo 

Moreover, the function v belongs to the class C 3 in some neighborhoods of the 
points x± and ±v'(x±) > 0. 

Note that if Assumption 12.21 is satisfied for some Ao , then it is also satisfied for 
all A from some neighborhood of Ao . 

Our goal in this section is to describe asymptotics as H — > of solutions u + (x) — 
u + (x;X,h) and it-(x) = U-(x;X, h) of equation (|2.1[) exponentially decaying as 
x — > +oo and x — oo, respectively. These asymptotics will be given in terms 
of an Airy function and are uniform with respect to x G [xi,oo) or x G (— oo,xi] 
where x\ is an arbitrary point from the interval (x_, x+). 

2.2. Let us recall the definition of Airy functions and their necessary properties 
(see, e.g., [7]), for details). Consider the equation 

-w"{t)+tw(t) =0 (2.2) 

and denote by Ai (t) its solution with asymptotics 

Ai(t) = 2- 1 7r- 1/2 i- 1/4 exp(-2i 3/2 /3)(l + 0(^ 3/2 )), H +oo. (2.3) 

Then 

Ai(t) -7r- 1 / 2 |ir 1 / 4 sin(2|i| 3 / 2 /3 + 7r/4) + 0(|tr 7 / 4 ), t ^ -oo. (2.4) 

Note that Ai (t) > for edit > 0. 

The solution Bi (t) of equation (|2.2p is defined by its asymptotics as i — >• — oo 
which differs from (|2.4|) only by the phase shift: 

Bi(i) = -^- 1 / 2 |tr 1 / 4 sin(2|i| 3 / 2 /3- 7 r/4) + 0(|ir 7/4 ), i -> -oo. (2.5) 

For f > 0, this function is positive and satisfies the estimate 

Bi (i) < C(l + t)" 1/4 exp(2t 3/2 /3). (2.6) 

Here and below we denote by C and c different positive constants whose precise 
values are of no importance. 

We also use that all asymptotics (|2.3[) . (|2.4p and (|2.5p can be differentiated in t. 
In particular, the Wronskian 

{Ai(t),Bi_(t)} := Ai'(t) Bi(t) - Ai(t)Bi'(i) = -tt" 1 . 
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It follows that 



Bi (s) Ai _1 (s) -Bi(t)Ai _1 (i) = tt" 1 ^ Ai " 2 (r)dr, s > t > 0. (2.7) 
2.3. To formulate results, we need the following auxiliary functions £,±{x) = 



/3 r x \ 2 / 3 



= -(2 \l(y)\ 1/2d y) > 



(2.8) 



and 



/3 \ 2 / 3 



Here is a list of properties of these functions. The following result is practically the 
same as Lemma 3.1 from Chapter 11 of [7]. 

Lemma 2.3. Let xi G TTiera £ + G C 3 (xi,oo), £_ e C 3 (— 00, xi) and 

£±(x) — > +00 as 1 -> ±00. TTie derivatives 

±&(x)>0, e±(x±)=±\v'(x ± )\^ 3 (2.9) 

and the functions £±(x) satisfy the equation 

&(x) 2 £±(z) =<?(*). (2-10) 

It follows from this lemma that the function 

P ±(x) = (i&(*)r 1/3 )"i&(*)r 3/a (2-n) 

is continuous. Moreover, using identity (|2.10[) . we see that 

- I6p ± (x) = ^ ± {x)- 2 +^ ± {x)^q{xr 2 q"{x) - 5q(x)- 3 q'(x) 2 ), x ± X±, (2.12) 
and hence according to Assumption 12.11 

\p±(x)\ < C\Z±(x)\-V a -", p = 3min{ Po -l,l}/2>0. (2.13) 

2.4. Let us construct solutions u±(x) = u±(x;\,h) of equation (12. ip with 
semiclassical asymptotics as h — > or (and) x — > ±00. We define these solutions 
by their asymptotics as x — > ±00. Below all asymptotic relations are supposed to 
be differentiable with respect to x. In this subsection, we only formulate results. 

Proposition 2.4. Under Assumvtion \2.1\ for every fixed h > 0, equation (|2.1[) has 
a (unique) solution u±(x) such that 

U±(x) =2~ 1 ir 1 / 2 h 1/6 q(x)- 1 / 4 exp ( =p tr 1 C q(y) 1/2 dy\ 
x (l + 0(\ f q{yf/ 2 dy\-^)) 

J x± 

where p\ = min{/9o — 1, 1} > as x — > ±00. 

Uniform asymptotic formulas for u±(x) are given in the following assertion. 
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Theorem 2.5. Let Assumptions \2.l\ and \2.2\ hold. Ifix > ±x±, then the solutions 
u±(x) = u±(x;\,ti) admit the representations 

u ± {x) = \e ± (x)\-^ 2 Ai(h- 2 / 3 U^)){l + e±(x;h)) (2.14) 

where the remainder satisfies the estimate 

\e±(x;\,h)\ <Ch{l + \Z ± (x)\)- p , p = 3min{p - l,l}/2 > 0. (2.15) 

Let Xi € (x—,x+). On the interval [xi,x+] (on the interval [x_,xi]) the function 
u+ (the function u_) admits the representation 

u±(x) = \e±(x)\- 1/2 Ai(h- 2 / 3 ^ ± (x)) + 0(h 7 ' & (h 2 ^ + \x - Xil)- 1 / 4 ). (2.16) 

Away from the points x±, we can replace the Airy function Ai (t) by its asymp- 
totics (|2.3j) or ()2.4j) . Indeed, in view of (|2.9|) . we see that 

|£±(x)| >c|x-x±|, c>0, (2.17) 

and hence ft~ 2 / 3 £ + (x) — > ±oo if hr 2 / 3 (x — x + ) — > ±oo and Fr 2 /' 3 ^- (x) — > ±oo if 
ft~ 2 / 3 (x — x_) — > =Foo. This leads to the following result. 

Corollary 2.6. Suppose that 5nh~ 2 / 3 > c > (in particular, 8% may be fixed). 
Then the functions u± (x) have asymptotics 

u ± (x) = 2- 1 ir 1 ' 2 h 1 '\(xr 1 ' i exp ( T h- 1 f X q(y) 1 ' 2 dy) (l + 0(h\£ ± (x)\- 3 / 2 )) 

(2.18) 

as h — > uniformly in x > x + + o/ /or n_|_(x) and in x < x_ — o/ /or u_(x). Let 
Xi G (x_,x + ). T/ien the functions u±(x) have asymptotics 

u±(x) = ^^^(x)!" 1 / 4 sin ( ± Tr 1 / X± |«(y)| 1/2 dy + tt/4) 

+ 0(fi 7/6 |x-x±r 7/4 ) (2.19) 

as ft — > uniformly in x £ [xi , x + — <5/J /or u + (x) and uniformly in x g [x_ + o^, Xi] 
/or w_(x). 

On the other hand, using estimates (|2.17|) and | Ai (t)\ < C(l + |t|) -1 / 4 , we obtain 
uniform in h estimates of the functions u±(x) in neighborhoods of the turning 
points. 

Corollary 2.7. For sufficiently small x — x±\, the estimate 

\u±(x)\ < C(l + h- 2 / 3 \x - xil)- 1 / 4 (2.20) 
holds with a constant C which does not depend on H. 

We note that all asymptotic relations (j2~T4"|) . (|2~T 7 J|) . (j2~T%|) and ((2~T9l can be 

differentiated with respect to x. In particular, we have asymptotics 

u'±(x) = Tvr 1/2 ft- 5/6 |g(x)| 1/4 cos ( ± ft" 1 \q(y)\ 1/2 dy + tt/4) 

+ 0(ft 1/6 |x-x±r 7/4 ) (2.21) 
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as ft — > uniformly in x £ [x\ , x+ — Sn] for u+ (x) and uniformly in a; £ [x- + Sn, x{\ 
for u-(x). All these relations can also be differentiated with respect to A. For 
example, we have 

du ± (x;\,h)/d\ = ±2- 1 Tr 1/2 h- 5/6 \q(x;\)\- 1/4 f ± \q(y; X)\- 1/2 dy 



X cos 



( ± h- 1 [ X± \q(y; X)\^ 2 dy + tt/a) + O^x - x ± \- 7 ^) 



as ft — > uniformly in x £ [x\ , x+ — <5/J for it + (a;) and uniformly in x £ [x- + Sn, X\\ 
for U-(x). 

2.5. Let us now calculate the norm of the function u± (x) in the space L 2 (x,\ , ±oo) 
where X\ £ (x_, x + ). Actually, we will obtain a more general result. 

Proposition 2.8. Let a function w(x) be differentiable on the interval (x^,x + ) 
except a finite number of points. Suppose that w{x)and w'(x) are locally bounded 
functions and that, for some N, 

[ x , N 

\w(x)\<Cq(x) / \q{y)\ 1/2 dy (2.22) 



if\x\ is large. Then under Assumvtions \2.\\ and \2.2\ we have the asymptotic relation 

rioo r'£± 

}{x)u 2 ± {x)dx = 2- 1 irh 1/3 / w(x)(X-v(x)y 1/2 dx + 0(h 2/3 ). (2.23) 

J X± 

Proof. We will prove (|2.23p for the sign " + " , omit this index and add h. Using 
asymptotics (|2.14D , (12. 15ft and the estimate > c > 0, we see that 

^:c+ + i poo 

w{x)u 2 h {x)dx < C £'(x)Ai 2 (h- 2 / 3 £,(x))dx = C 1 h 2 / 3 . 
f •' x + 
Similarly, using identity (|2 . 10[) and condition (12 .22[) we find that 

/>oo 

,{x)u 2 h {x)dx <C Z'(x)£ 3N / 2+1 (x)Ai 2 (h- 2 / 3 Z(x))dx = 0(h°°). 

Ix + + l Jx++1 

Suppose that 5n — > as h — > but Shhr 2 ^ 3 > c > 0. The integral of u\{x) 
over {x+ — 5fr,x + ) is estimated by C5% because according to (|2.20[) the functions 
un(x) are uniformly bounded in a neighborhood of the point x+. On the interval 
(xi,x + — Sn), we have a relation 

i-x + —6 n rx + -S h 

/ w(x)u 2 l (x)dx = irh 1 ^ 3 / w(x)\q(x)\ 

J x-[ J x-i 



1-1/2 



(ft- 1 f + \q(y)\ 1/2 dy + n/A)dx + 0{h i ' 3 5 h 1 ). (2.24) 



x sin 2 ' 

Indeed, in view of asymptotics (|2.19p we have to show that the integrals 

ft 7 / 3 /" + ^ \x~x + \- 7 ' 2 dx and ft 4 / 3 5h ^(x^'^x-x+^^dx 

J X\ J X\ 

arc 0(h 4 / 3 S h x ). The first of them equals Ch 7 / 3 S^ 5/2 which is C^ft 4 / 3 ^ 1 ) because 

ft = 0(Sf/ 2 ). To estimate the second integral, we have to additionally take into 
account that 

\q(x)\>c(x+-x), c> 0. (2.25) 
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Next, we replace sin 2 (-) in the right-hand side of (|2.24|) by 1/2. Let us esti- 
mate the error. Integrating by parts separately on every interval where w(x) is 
differentiable, we see that 



* w(:r)|<7(a:)|~ 1/2 exp^i/T 1 j + \q{y)\ l ' 2 dy^dir 
2- 1 ihw(x+ - S h )q{x+ - Sh)- 1 exp^TT 1 f " \q(y)\ 1/2 dy) 



x+Sn 



r x + 

I x + -8 h 

+2~ 1 ih I (w'{x)q(x)- 1 - v'(x)q{x)- 2 w(x)) 

J Xi 

xexp^ih- 1 J + \q(y)\ 1/2 dy^dx + 0(h). 

The right-hand side here is bounded by 

Ch(l + \q(x+ - Sh)]- 1 + jf + K q(x)- 2 dx) 

which in view of estimate ()2.25|) does not exceed CHS^ 1 . Thus, it follows from 
(pQ4l) that 

w(x)u 2 h (x)dx = 2- 1 nh 1/3 [ + * w(x)(X - v(x)y 1/2 dx + O^^d^). 



Finally, making an error of order 0(h 1 ^ 3 6^ 2 ), we can extend the integral in the 
right-hand side to the whole interval (x\ : x + ). Setting 5n = h 2 / 3, and putting the 
results obtained together, we arrive at asymptotic relation (|2.23[) . □ 

Of course (|2.22|) is a very mild restriction. It is satisfied for v(x) — w(x). It is 
also true for all functions v(x) if w(x) if bounded by some power of \x\ at infinity 
and is even less restrictive if v(x) — > oo as \x\ — > oo. In particular, setting w(x) = 1, 
we obtain 

Corollary 2.9. The asymptotic relation holds: 

f±00 fX± 

u 2 ± {x)dx = 2- 1 7rn 1 / 3 / (X-v(x))-^ 2 dx + 0{h 2 ^). (2.26) 

I X\ J X\ 

3. Proof of Theorem [231 

3.1. We will prove Theorem 12.51 for the sign " + " and omit this index. On the 
contrary, we add the index h to emphasize the dependence on it of various objects. 
Let X\ £ (x_,x+), x £ (xi,oo) and let the function £(x) be defined by formulas 
(|2.8[) . According to Lemma l2~3l £ (x) £ (£i,oo) where £i = £(xi) and x can be 
considered as a function of £ if £ £ oo). 

Let us make the change of variables x n> £ in equation (|2.1j) and set 

u h {x) = e(x)- 1/2 Mh- 2/3 a*))- (3-1) 
Then using identity (|2.10l) . we obtain that 

- mn- 2/3 + h- 2 /^Mh- 2 / 3 o = h^r(Oh(h- 2 / 3 o, (3-2) 

where 

r(0= P (x{0) (3.3) 
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and p(x) is defined by formula (|2.1iJ . In view of (|2.13|) we have the estimate 

KOI <C(l + |eir 1/2 - p , p = 3min{p -l,l}/2>0. (3.4) 

Setting in p.2[) t — ft~ 2 / 3 £, we get the following intermediary result. 

Lemma 3.1. Let t = h~ 2 ^ 3 ^(x), and let the functions uti(x) and fn(t) be related 
by formula (|3.ip ■ Then equation (|2.ip for x > x\ is equivalent to the equation 

-fn(t)+tf h (t) = R h (t)f h (t) for t>i x tr 2 ' 3 (3.5) 

where 

R H (t) = ft 4 / 3 r(ft 2 / 3 i). (3.6) 

3.2. Let us reduce differential equation (|3.5I) to a Volterra integral equation. Set 
K h {t,s) = -TT(Ai(t)Bi(s)-Ai(s)Bi(t))R h (s), s>t, (3.7) 
and consider the equation 

f h {t) = Ai (t) + / K h (t, S )f h (s)ds. (3.8) 



Differentiating it twice, we see that its solution satisfies also differential equation 
(|3.5p . We will study equations (|3 . 5|) or (|3.8[) separately for t > and t < 0. 



Lemma 3.2. For t > 0, equation (|3.5p /ias a solution fn(t) such that 

f h (t) =Ai(t) (l + %(i)) (3.9) 

w/iere 

|%(t)| < C^(l + ^ 2/3 t)-", p = 3min{p - l,l}/2 > 0. (3.10) 
Proof. Making the multiplicative change of variables 

/»(*)= Ai(t) 5B (t) (3.11) 
and using (|2.7[) . we rewrite equation p.8p as 



3ft (t) = l- y L n (t,s)g A (s)ds, (3.12) 

where 

L n (t,s)=Ai(t)- 1 K H (t,s)Ai{s)= J Ai ~ 2 (r)dT Ai 2 (s)R h (s) , s > t. 
It follows from (|2~3|) that 

J Ai " 2 (r)dr < C*exp(4s 3 / 2 /3) 



so that according to (|3 .4[) and (|3.6I) 

|L fc (t, «)| < Cn 4/3 s- 1/2 (l + h 2/3 s)- 1/2 - p , 0<t<s. (3.13) 

This estimate allows us to solve equation (I3.12j) by iterations. In particular, the 
solution of (|3.12p satisfies the estimate 



\gn(t)-l\<Cj \L h (t,s)\ds. 

Now estimate (|3.10p on the remainder rjn(t) = gnif) — 1 follows again from (|3.13p . 

□ 
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Putting together formulas p.lj) and (|3.9|) . we obtain representation (|2.14j) with 
e h (x) = rjfi(h~ 2 / 3 !;(x)). Estimate (|3.10[) implies estimate (|2 . 15[) . This leads to 
the assertion of Theorem 12.51 for x > x + . In particular, for a fixed h, we get 
Proposition ^. 41 

Next, we consider the case t < 0. 

Lemma 3.3. For t G [fi,~ 2 / 3 £i, 0], i/ie solution fn(t) of equation (|3.5|) satisfies the 
estimate 

\Mt)-Ai(t)\<Ch(l + \t\)-^ 4 . (3.14) 
Proof. Let us rewrite equation (|3.8[) as 

/»(*) = /fi 0) («) + f K h (t,s)f H (s)ds, (3.15) 



'a 

where the new "free" term 

/>oo 

/< 0) (t) = Ai (t) + /«(*) = / Jjf ft (t, S )/ ft ( S )ds. (3.16) 

Jo 

It follows from JOJ and (JUJ) that 

Ai 2 (t)+Ai(i)Bi(i) < C*(l + ^ 1/2 , t > 0, 
and from (j2~4|) and (j2~B1 that 

|Ai(t)| + |Bi(t)| <C(l + |i|)" 1/4 , i<0. (3.17) 
Therefore using f|3T4|) . (|3T6]) and (|3~7]) . we find that 

/"OG /»oo 

\f£\t)\ < C(| Ai (t)| jf Ai(«)Bi(«)|iJfc(«)|ds + |Bi(t)| jf Ai 2 (s)|i? ft ( S )|eZs 

< Cl(l + Itlr 1 / 4 ^/ 3 / S -V2(l + ft 2 /^)- 1 / 2 -^ 

Jo 

<C 2 (l + \t\)- 1/4 h. (3.18) 

Let us now consider equation ()3.15p . By virtue of estimates (|3.4j) and (|3.17p its 
kernel satisfies the bound 

\K h (t,s)\ <Cft 4/3 (l + |t|)- 1/4 (l + |s|)- 1 / 4 r(n 2 / 3 s), i<s<0, (3.19) 

where the function r(fi, 2 / 3 s) can be estimated by a constant. Thus, solving (j3. 15[) 
again by iterations, we obtain the estimate 

\fh(t) - fi°\t)\ < C x J° \K h (t,s)\(l + \s\)- l / 4 ds < C 2 h^ 3 f(l + | s |r 3/4 d* 

< C 3 ^ 4/3 (l + |t|) 1/4 . (3.20) 

If t 6 [fi~ 2/3 £i, 0], then combining definition (f3~16l) with estimates ([3~18| and (|3~20| . 
we get estimate p. 141) . □ 



In view of formula (|3.1|) , this lemma yields the result of Theorem 12.51 for x G 

Differentiating integral equation (|3.8|) with respect to i, we obtain asymptotic 
relations for f' h (t) and then for u' h (x). This concludes the proof of Theorem 12.51 
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4. Semiclassical ASYMPTOTICS OF EIGENFUNCTIONS 

4.1. Let A = A (ft) be an eigenvalue of the Schrodinger operator Hn = —h 2 d 2 /dx 2 + 
v{x) from a neighborhood of a non-critical point Ao- Then the solutions u±(x) are 
proportional: 

u_(x; A, ft) — a(A, h)u + (x; A, ft). (4-1) 

Choose an arbitrary interior point x of the interval (x_(A),a;+(A)). To calculate 
the Wronskian of u+{x) and U-(x), we use asymptotic relations (|2.19p and (|2.21l) . 
Setting 

<p ± (x;\) = ± (X-viy^dy, x G (a;_(A),a: + (A)), (4.2) 

J X 

we find that 

w(X, ft) — u + (x] A, h)u'_(x; A, ft) — u-(x; A, h)u' + (x; A, ft) 

= 7rft~ 2/3 (sin(ft~V+0c; + t/4) cos(ft~ V- (»; A) + 7r/4) 
+ cos(ft _ V+(aj; A ) + "V 4 ) sinih^ip-ix; A + tt/4)) + <3(ft 1/3 ) 

= irh- 2 ' 3 sin(^ 1 $(A) + tt/2) + ©(ft 1 / 3 ) (4.3) 
where $(A) = (p+(x; A) + <y3_(x; A) so that 

rx+{\) 

$(A) = / (A - v{y))^ 2 dy. (4.4) 

J X-{\) 

Since u>(A, ft) = 0, we see that 

sin(ft- 1 $(A) +71-/2) = 0(h) 

and hence 

l-X+{\) 

/ (A-v(ir)) 1/2 d:r = 7r(ri + l/2)ft+0(ft 2 ) (4.5) 

for some integer number n — n(A, ft). This gives us the famous Bohr-Sommerfeld 
quantization condition. 

Suppose now that a number ir(n + l/2)ft belongs to a neighborhood of Ao- Let 
us check that there exists an eigenvalue A„(ft) of the operator Hn satisfying the 
estimate 

|$(A„(ft)) -7r(n+ l/2)ft| < Cft 2 . (4.6) 
Since u± G L 2 (R±), it suffices to show that w(X, ft) = for some A = A„(ft) 
satisfying estimate (|4.6p . Using the equality A — v(x±(X)) — 0, we find that 

rx + (\) 

$'(A) = 2- 1 / (X-v(y)y^ 2 dy>0. (4.7) 

Hence $ is a one-to-one mapping of a neighborhood of Ao on a neighborhood of 
= ^(Ao)- Set [i = $(A) and 

e(/i, ft) = vr" 1 ft 2/3 w($" 1 (^), ft) - sin(ft"V + i"/2)- (4.8) 

In view of (|4.3I) this function satisfies the estimate |e(^t, ft)| < Cft with a constant 
C which does not depend on ft and fi from a neighborhood of /j,o- We have to show 
that the equation 

sin(ft~ V + tt/2) + e(/i, ft) = 
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has a solution /i n (ft) obeying the estimate 

\fl, n (h) -7r(n + l/2)ft| < Ch 2 . 

Setting s = hr 1 /i + tt/2, we see that this assertion is equivalent to the existence of 
a solution s = s n (h) of the equation 

sin s + e(ft(s - tt/2), ft) = (4.9) 

obeying the estimate 

|s n (ft)-7r(n + l)| <Cft. (4.10) 

The last fact is obvious because e(h(s — 7r/2), ft) = O(ft). 

Next, we will show that for every n there is only one eigenvalue of the operator 
Hfi satisfying (|4.6I) . To that end, we have to check that equation (j4.9[) cannot have 
two solutions satisfying (|4. 10[) . Supposing the contrary, we find a point s = s n (h) 
such that 

de 

coss = -ft— (ft(s-7r/2),ft) (4.11) 

and s n (h) = ir(n + 1) + 0(K). Observe that relation (|4.3[) can be differentiated in 
A which yields 

dw(X, h)/d\ = vr$'(A)ft- 5/3 cos(ft _1 $(A) + tt/2) + 0(?T 2/3 ). 

It follows that function (|4.8|) obeys the estimate de(/i,H)/dfi — 0(1). Thus, the 
right-hand side of equation (j4.11j) is 0(h) while its left-hand side tends to (— 1)™ +1 
as h — >• 0. 

Finally, plugging asymptotics (|2 . 19[) and (I2.2ip into (|4.ip . we find that 

sin(ft-V-(»; A) + tt/4) + O(h) 

=a{\,h)(sm(h- 1 ip+{x;\) + vr/4) + 0(h)) (4.12) 

and 

cos(7rV-0; A ) + V 4 ) + O(ft) 

= - a(A, ft) ( cos(ft- V+(x; A) + vr/4) + O(ft)) . (4.13) 

Together, these two relations imply that |a(A, ft)| = 1 + O(ft). Moreover, since 

ip+(x; A) + ^_(a;; A) = n(n + l/2)ft + 0(ft 2 ), 

it follows from (gTTgl and P~llfl) that 

a(A,ft) = (-l)" + 0(ft). (4.14) 

Thus, we have obtained the following result. 

Theorem 4.1. Let Assumptions and 12.21 ftoM /or a point Xq. Suppose that 
an eigenvalue A = A(ft) of the operator belongs to a neighborhood of Xq. Then 
necessarily condition (|4.5[) is satisfied with some integer number n — n(X, ft). Con- 
versely, for every n such that 7r(n + l/2)ft belongs to a neighborhood o/$(Ao), there 
exists an eigenvalue X n (h) of the operator Hn satisfying estimate (|4.6j) with a con- 
stant C not depending on n and ft. Such an eigenvalue A„(ft) is unique. Moreover, 
the coefficient a(A, ft) in (I4.ip has asymptotics (|4.14[) where n is the same number 
as in (|4.5[) . 
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Corollary 4.2. Let an interval (01,02) belong to a neighborhood of a point Ao 
satisfying Assumptions [2~T1 and 12.21 Then the total number Nn of eigenvalues of 
the operator Hji in this interval equals 

N h = 7r- 1 ($(o 2 ) - ^{ ai ))h- 1 + e(H) (4.15) 

where |e(7l)| < 1 for sufficiently small h. 

Proof. According to Theorem 14.11 there is exactly one eigenvalue of the operator Hn 
in a neighborhood of size Ch 2 of every point <& _1 (7r(ra+l/2)?i). These neighborhoods 
have empty intersections for sufficiently small h. Thus, iV^ equals the number of 
points ir(n+ 1/2)K lying in the interval ($(ai), < &(a2)). Clearly, this number equals 
the right-hand side of (|4.15|) . □ 

Remark 4.3. Suppose that Assumptions 12.11 and 12.21 hold true for all A G [01,02]. 
Then remainders in different asymptotic formulas of this paper can be estimated 
uniformly in A € [01,02]. Formula (|4.15[) also remains true for such (01,02). 

Note that definition (|4.4[) can be rewritten as 

$(A) = T 1 / / dpdx. (4.16) 

J Jp 2 +v(x)<\ 

Indeed, integrating in the right-hand side first over p we obtain the right-hand side 
of (|4.4[) . It follows that the asymptotic coefficient in (|4.15l) is the volume of a part 
of the phase space: 

<J>(o 2 ) ~ $(oi) = 2 _1 mes{(x,p) <E R 2 : di < p 2 + v(x) < o 2 }. 

Thus, relation (|4.15j) is the semiclassical Weyl formula with a strong estimate of 
the remainder. 

4.2. Let us denote by ijj{x) — tp(x; A, H) the eigenfunction of the operator Hn 
corresponding to its eigenvalue A. We suppose that ifi — i\> £ L 2 (R + ) and ||^|| = 1 
which fixes ip up to a sign. Clearly, 

ip(x) = c±u±(x), c±=c±(\,h), (4.17) 

where according to (|4.14p 

\c+(x,h)\ = \c-(\,h)\(i + o(h)). 

Therefore it follows from Corollary 12.91 that 

fX+(X) -1/2 

\c ± (\,h)\=2 1 ^n- 1 ^h- 1 ^( (A - v{x))~ x / 2 dx) +O(o 1 / 6 ), (4.18) 

V Jx_(A) ' 

which in view of Theorem 12.51 yields the following result. 

Theorem 4.4. Under the assumptions of Theorem 14.11 let us denote by tp(X, K) 
the real normalized eigenfunction (defined up to a sign) of the operator Hti corre- 
sponding to its eigenvalue A = X(h). Let X\ be an arbitrary point from the interval 
(a;_(A), x + (X)). Then, for x € (xi,oo), asymptotics of ip(x] A, K) as h — » is given 
by formulas (|4.17l) . (14. 18)) for the sign "+" and asymptotic relations of Theorem \l.h\ 
for the function u+(x; A, H). Similarly, for x € (—00, xi), asymptotics of ij)(x; A, ti) 
as h — > is given by formulas (|4.17[) . (14. 18)) for the sign " — " and asymptotic rela- 
tions of Theorem \2.5\ for the function U—(x;X,K). In neighborhoods of the turning 
points, the estimate 

\ip(x; A, h)\ < C(h 2 / 3 + \x- x±|)- 1/4 (4.19) 
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holds with a constant C which does not depend on h. 

In view of formula (I2.23|) . this result can be supplemented by the following 

Proposition 4.5. Let a function w satisfy the assumptions of Proposition 12.81 
Then under the assumptions of Theorem 14.41 we have 

w(x)ip 2 (x; A, h)dx = / w(x)(X — v{x))~ 1 / 2 dx 

x(/ (A - v{x))~ l l dx) +C>(ft 1/3 ). (4.20) 

Vj_(A) ' 

In particular, this relation applies to the potential energy 

/oo 
v(x)i> 2 (x;\,h)dx 
-oo 

and by virtue of the energy conservation K(X, K) + V(X, H) = A, we also obtain the 
asymptotics of the kinetic energy. 



Corollary 4.6. Under the assumptions of Theorem 14.41 asymptotic relation (|1.3p 
holds with the leading term K c i (A) given by (jl.4|) . 

Since K c i(X) > 0, for small H the kinetic energy K(X, H) > c > or, equivalently, 
the potential energy V(A, ft) < A — c. This implies that the eigenfunctions ip(x; A, h) 
are not too strongly localized in neighborhoods of the turning points x± (A) . In view 
of estimate (|4. 19[) . this statement can be reinforced. 

Proposition 4.7. Let the assumptions of Theorem \4.1\ hold. and let \\ip(\,h)\\ = 1. 
Then 

r x ± (X)+8 

/ i[> 2 (x;X,H)dx < CS 1/2 

where the constant C does not depend on H. 

Remark 4.8. It follows from asymptotics (|2.14l) and (|4.18[) that 

|V(ar±(A); A, ft)| = a±(\)h~ 1 / 6 (l + 0(^ 1/3 )), 
where the coefficient 

, / r x+< -^ , \-i/2 

a±(A) = 2 1/2 tt- 1/2 (J (X-v(x))- 1/2 dxj |t/(x±(A))r 1/6 Ai(0) ^ 0. 

This contradicts the assertion of Theorem 7.1 of [8] that normalized eigenfunctions 
are uniformly bounded in neighborhoods of turning points. 

4.3. Recall that a classical particle (of mass m and energy A) moves periodically 
(see, e.g., [6]) in a potential well bounded by the points x_ = x_(A) and x+ = a; + (A) 
such that v(x±) = X. Let us check that the asymptotic coefficient K c i in (|1.3|) 
coincides with the averaged over the period T — T(X) value 



K av = T- 1 f K{t)dt 
Jo 



of the classical kinetic energy 

K(t) = mx'{i) 2 /2 = A - v[x{t)). 
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Since 

dt = x'(t)- l dx = (m/2) 1/2 (A - v(x))- 1/2 dx, 
the period is given by the formula 

r x + dt , f x + 

T = 2 ^dx = (2m) 1 / 2 / (X-v(x))-^ 2 dx 

J X - dx J x _ 

and 

r T r x+ 

K av — T^ 1 / (X - v(x(t)))dt = 2(m/2) 1/2 T- 1 / (X - v(x)) 1/2 dx. 



Jo J x- 

Putting these two relations together, we obtain for K av the same expression (II .4[) 
as for K c \. This proves the equality 

K cl (X) =K av (X). 

We also note that 

K cl {\) = (2dln$(A)/dA) _1 
where the function $(A) is defined by formulas (I4.4[) or, cquivalently, (|4.16[) . For 
the proof of this equality, it suffices to plug representation (|4.7|) for the function 
$'(A) into formula (fOl). 

5. Discontinuous potentials 

5.1. Away from the turning points, assumptions on v(x) can be somewhat 
relaxed. Consider, for example, an interval (x- + S,x+ — S) where S > 0. There, it 
suffices to require that v G C 1 and that v' be absolutely continuous so that v" G L 1 
(instead of v £ C 2 ). In this case the function r(£) defined by formulas (|2.11[) and 
(|3.3j) belongs to L 1 only so that the factor r(fi 2 / 3 s) in the right-hand side of estimate 
(13.191) cannot be neglected. Therefore (cf. (|3.20[1 ) we have the estimate 

\K H (t, s)|(l + Isl)- 1 / 4 ^ < Ch 4 ^(l + Itl)- 1 / 4 / \r(h 2 / 3 s)\ds 

<Cin 2/3 (l + |<|)- 1/4 / \r(s)\ds. 

It follows that instead of p,14p we have a slightly weaker estimate with fi 2 / 3 in 
place of h in the right-hand side. All other estimates remain unchanged. Thus, 
Theorem 12.51 is true with a little bit weaker estimates of remainders in asymptotic 
formulas for u±(x] X, h) inside the interval (x_+<5, x+ — S). Repeating the arguments 
of Section 4, we get the following result. 

Proposition 5.1. Under the assumptions above, all results of Theorem lA.lK and of 
Corollarv \4,.2\ about eigenvalues of the operators Hn remain true with the remain- 
ders 0(h 5 / 3 ) in ((431), 0(ft 2 / 3 ) in (|4~T4|) and Ch 5 / 3 in the right-hand side of (1431 . 
Theorem \AA\ about corresponding eigenf unctions remains also true with the remain- 
ders 0(h 5 / 6 \x - x±\- 7 / i ) in (gZEaJ) , 0(h- 1 / 6 \x-x±\- 7 / i ) in (l2~2lT) and Oih 1 / 18 ) 
in (|4TTgj) . 

5.2. Our goal in this subsection is to extend the results of Section 4 to functions 
v(x) with a singular point xq inside a potential well. 

We suppose that Assumption 12.11 holds everywhere except a point xq and that 
Assumption 12.21 holds for some Aq such that Xq is an interior point of the interval 
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(x-(Xq), x+(Ao)). We assume that v[x) has finite limits at xo but the left and 
right limits might be different. Finally, we require that v' € L 2 (xq,xo ± S) and 
v" E L 1 (xo,x ± 5) for some 5 > 0. 

Now we can construct solutions u + (x) and u_ (a;) of equation (jl.ip on the in- 
tervals (xo,oo) and (— oo,xo), respectively. Define, as usual, the function r(£) by 
formulas (|2.11[) and (j33J). Since r £ L 1 (xo,xo ± 5), the limits w±(a;o ± 0) and 
u±(xq ± 0) exist, and we can use formulas (|2.19|) and (|2.21j) for these limits (with 
slightly weaker estimates of the remainders - see subs. 5.1). It follows that the 
Wronskian w(X, H) of w + and u_ calculated at the point xq is given by the expres- 
sion (cf. flOJ)) 



Let an eigenvalue A of the operator H h be close to A . Since w(X, h) — 0, we see 
that 

p{xq, X) sin(ft _1 ^ + (xo; A) + 7r/4) cos{h~ x ip-{xa; A) + 7r/4) 
+p(x ,A)- 1 sin(?rV+OEo;A) +7r/4)cos(n - V-(a:o;A) +tt/4) = 0(fr 2 / 3 ). (5.2) 

Formula (15.2[) yields a generalization of the Bohr-Sommerfeld quantization condi- 
tion (|4.5|) and reduces to it if u(xo + 0) = v(xa — 0). 

Let the coefficient a(A, be defined by equality (|4.ip . To calculate |a(A, fi)|, we 
use again relations (14. 12|) and (|4. 13[) . However, additional factors \q(xo — 0)| -1 / 4 
and \q{xo — 0)! 1 / 4 appear now in their left-hand sides. Similarly, additional factors 
\q(x + 0)|~ 1 / 4 and \q(x + 0)! 1 / 4 appear in their right-hand sides. This implies that 



a 2 {X,K) =p 2 (x ,X)cos 2 (h V-( s o; A) + 7r/4) 

+ p- 2 (x , A) sin 2 (ft-V-(^o; A) + tt/4) + 0(h 2 ^) 
= (p 2 (x ,X) sin 2 {IT V+ [xo ; A) + tt/4) 

+ p" 2 ^, A) cos 2 (n-V+(^o; A) + tt/4))" 1 + 0{h 2/3 ). (5.3) 



As before, using formula (|2.26l) (where x\ = xo) and the normalization condition 
| 1 1 = 1, we obtain explicit expressions for the absolute values of constants c±(A, H) 




where 




(5.1) 



in (14171) : 




(5.4) 



and 




(5.5) 
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Thus, Theorems 14.11 and 14.41 can be supplemented by the following result. 

Theorem 5.2. Under the assumptions above, let an eigenvalue A = X(h) of the 
operator Hn belong to a neighborhood of Xq. Then necessarily condition (|5.2p is 
satisfied with the numbers f±(xo;X) and p(xo, X) defined by (|4.2[) and (|5.1[) . re- 
spectively. All assertions (for x\ — Xq) of Theorem 14.41 about the corresponding 
normalized eigenfunction ip(x;X,h) are true with the constants c±(X,h) whose ab- 
solute values are determined by formulas (|5.3[) . (|5.4|) and (|5.5|) . 



Remark 5.3. If the functions v'(x) and v"(x) are bounded in a neighborhood of 
the point xq, then even in the case v(xq + 0) ^ v(xq — 0) estimates of all remainders 
are the same as in Section 4. Thus, we have 0(h) in (|5.2p and (^(/i 1 / 6 ) in ()5.3[) - 
(ESI). 



Remark 5.4. Let under the assumptions above v(xq + 0) = v(xq — 0). Then 
all conclusions of Proposition 15.11 remain true although the function v'(x) is not 
required to be continuous at the point xq- In particular, we see that jumps of 
derivatives of the function v(x) at the point xq are inessential. 

5.3. Let us consider an explicit example: 

v(x) — a + + v + x a+ for x > and v(x) = a_ + w_ \x\ a ~ for x < 0, (5.6) 

where v± > and a± > 0. Then all A > max{a + , a_} are non-critical, the equation 
v(x) — X has two solutions x+ > 0, x_ < and (x-,x+) is a potential well. The 
point xq — might be singular and p(0, X) = (X — a-) 1 / 4 (A — a + ) -1 / 4 . 

For potentials (|5.6|) . the integrals in formulas (|4.4|) and f)4. 18[> can be calculated 
in terms of the beta function B. Observe that x + = (Ad -1 ) 1 /" if v(x) = vx a for 
x > 0. For the integrals over (0, x+), we have 

/^(A - vx a y/ 2 dx = X 1/2 (X/v) 1 / a a- 1 B(3/2, 1/a) (5.7) 
Jo 

and 

/^(A - vx a )-^ 2 dx = A- 1/2 (A/w) 1/Q a -1 B(l/2, 1/a). (5.8) 
Jo 

The integrals over (x-,0) can be calculated quite similarly. 
It follows that the quantization condition (|5.2[) holds with 

9+ (0, X) = Ay 2+1/ « + u ; 1/a+ a; 1 B(3/2, l/a+), 

<^_(0, A) = A 1/2+1/Q -C 1/a -aI 1 B(3/2, l/a_), 

where A± = A — a±. In particular, in the case a + = a_ =: a the Bohr-Sommerfeld 
quantization condition reads as 

(A - a) 1 / 2 + 1 /"+i;; 1/a+ a; 1 B(3/2, l/a+) 

+ (X- a) 1 / 2+1 ^ a -vZ 1/a -aZ 1 B(3/2, 1/a-) = nh(n +1/2) + 0(h 2 ). (5.9) 

Plugging expressions (|5.7p and (|5.8p into (|1.4p we also find that 

X 1 / 2 + 1 / a + v - a + a -}B(3/2, l/a+) + X 1 J 2+1/a -vZ a -aZ 1 B(3/2 1 1/a.) 



K cl (X) 



A7 1/2+1/a+ «7 Q+ a7 1 B(l/2, l/a + ) + AI 1/2+1/ "-i;: a -a: 1 B(l/2, l/a_) ' 



Observe that Theorems 14. II and 14.41 can be applied to potential (|5.6|) if a+ = a_ 
and a± > 2. If a+ = a_ but a± £ [1,2), then we have to use Remark 15 A\ in 
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this case 0{h 2 ) in jOJ should be replaced by 0{h 5 / 3 ). If a± are arbitrary and 
a± > 1, then the conditions of Theorem 15.21 are satisfied. Moreover, according to 
Remark 15.31 in the case a± > 2, the estimates of the remainders can be improved. 
Finally, we note that if ay < 1, then v" $ ^(—5,6) so that the semiclassical 
approximation does not directly work (even for a+ = a_) although all formulas 
above remain meaningful. 

5.4. Let us briefly consider the problem on the half-axis. We now suppose that 
equation is satisfied for x > 0, tp € £ 2 (K+) and V'(O) = 0. Assumptions 12. ll 
and [22] should be slightly modified. Namely, we assume that the equation v(x) = A 
has only one solution x+ = x+(A) and v'(x+) > so that (0, x+) is a potential well. 
We suppose that the limit of v(x) as x — ► exists and that the functions v'(x) and 
v"(x) are bounded in a neighborhood of x = 0. Then the results of Theorem l2.5l on 
the solution u+(x) of equation (|1.1[) are true for all x > 0. In particular, it follows 
from formula (|2.19[) that 

rx+(\) 

u+(0;X,h) = TT 1 / 2 h 1 / 6 (\-v(0)y 1/i sm(h- 1 / (\-v(x))V 2 dx+w/4)+0(h r / 6 ). 

Jo 

(5.10) 

Since ip(x; A, = c+m+(x; A, K), this yields the quantization condition 

x+(A) 

(A - w(x)) 1/2 dx = rr^n + 3/4) + 0(h 2 ) 

where n = n(A, K) is an integer. 

Consider now the boundary condition ^'(0) = b^(0), b — b. It follows from 
(pT2Tj) that 

u' + (0] A, ft) = -tt^^-^^A - v(0)) 1/4 cos (h- 1 / (A - ^(x)) 1 / 2 ^ + rr/4) 

Jo 

+ o(h 1/6 ). 

Comparing this formula with (|5 . 10[) . we see that the value of w+(0; A, H) is inessential 
so that the quantization condition looks like 

x+(X) 

(A - v(x)) 1/2 dx = 7rh(n + 1/4) + 0(h 2 ). 

It does not depend on b. 

Other results of Section 4 can also be naturally extended to the problem on the 
half- axis. 

I thank D. Robert for a discussion of papers on microlocal analysis. 
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